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A direct method is proposed for reconstruction of the effective atomic 
number by means of multi-energy radiography of the material. The accuracy 
of the method is up to 95% . Advantages over conventional radiographic 
methods, which ensure accuracy of just about 50% , are discussed. A physi- 
cal model has been constructed, and general expressions have been obtained 
for description of the effective atomic number in a two-energy monitoring 
scheme. A universal dependence has been predicted for the effective atomic 
number as a function of relative (two-energy) radiographic reflex. The es- 
tablished theoretical law is confirmed by the experimental data presented. 
The proposed development can find multiple applications in non-destructive 
testing and related fields, including those in the civil sphere as well as anti- 
terrorist activities. 



07.85.-m ,78.20.Bh ,81.70.-q, 81.70.Fy ,87.59.Hp 



1 



1. Introduction 



Non-destructive testing (NDT) of different unknown materials, objects 
and media, including complex micro- and macro-structures, is an important 
field in modern applied physics 1J. Physical methods of NDT find multiple 
applications in a wide variety of technical areas, including quality control of 
industrial products, inspection of welded joints and connections, and con- 
firmation of building construction, pipelines, etc. An important role is also 
played by NDT in radioactive materials control |2] , security protection in avi- 
ation [3], jl], railway and automobile transport, and customs inspection of 
goods and luggage. In addition, broad applications related to the use of new 
NDT technologies in medicine [S], j^j, jZjj [Eh including separate diagnostics 
of tissues and organs, can be expected. Among the promising directions of 
NDT, one should note digital (computer) radiography. In it, digital recon- 
struction of an image of an object examined by X-ray or gamma-radiation 
is prepared based on the data collected in real-time detection. The most 
widely used radiation detectors for such applications are solid-state ID and 
2D scintillation systems 0, [10], [II], [12], [12] 

The main task of digital radiography, as well as of other NDT methods, 
is reconstruction of the physical parameters that determine technical char- 
acteristics and properties of the controlled objects being examined. This is 
traditionally carried out by reconstruction of the spatial configuration, ob- 
serving small details and defects of the image, and determining constitutive 
components and dimensions (thickness). In recent times, strong interest has 
been emerging in the development of direct methods for reconstruction of 
the overall composition of the object, i.e., quantitative determination of the 
effective atomic number, density, chemical composition, etc. To a large ex- 
tent, this is due to improvements in technical capabilities for detection and 
image processing and is driven by the desire for additional detail in the im- 
age resulting from the collected data. Often it is not sufficient to know just 
the geometrical structure of an object to ensure its complete identification 
for decision-making purposes. For example, fully reliable detection of ex- 
plosives is not possible without discrimination of materials by their effective 
atomic number (organics and inorganics), density, content of predominant 
chemical elements, etc. The importance of such enhanced monitoring cannot 
be overestimated. It is vitally necessary not only to provide enhanced tools 
for scientific and technological investigations, but to meet current needs for 
improved protection against terrorist threats to the safety and health of civil 
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populations. An important milestone on the way to solving this problem is 
the implementation of two-, and in general, of multi-energy radiography [T3] . 

m- m- en- 

2. Effective atomic number 

Among the parameters determining the constitutive structure of an un- 
known object or material, one should especially note the effective atomic 
number Z e s. In fact, this measure can provide an initial estimation of the 
chemical composition of the material. A large Z e g generally corresponds to 
inorganic compounds and metals, while a small Z c g < 10 is an indicator 
of organic substances. For many applications, e.g., radioisotope monitoring, 
cross-section studies of absorption, scattering and attenuation of electromag- 
netic radiation, testing of multi-component, heterogeneous and composite 
materials, etc., this parameter is of principal significance. It is very impor- 
tant to obtain its quantitative evaluation from the data provided by radio- 
graphic measurements. With the aim of making this possible, in the present 
work we have developed a physical model for two-energy radiography as a 
particular case of the multi-energy radiography. We show that a multiplicity 
of M = 2 is sufficient for direct determination of Z e g. General theoretical 
expressions are provided for determining the effective atomic number with 
respect to the photo and pair generation effects. 

The method proposed here for determination of Z e g is substantially dif- 
ferent from the approaches used previously. In these procedures, the effective 
atomic number is comprised of a linear superposition of several reference val- 
ues, i.e., it is expressed in the form of Z e g = J2k=iPkZk, where J2kPk = 1, 
and iV = 2,3,... corresponds to the multiplicity of the radiography used. 
Such "synthesis" of Z e g is quite similar to the formation of light of arbitrary 
color from several basic colors (Fig.l). Weights pk, or the relative content of 
the components (light Z^, he&vyZn, intermediate-Z^f), are determined from 
the data of radiographic measurements. Theoretically, this corresponds to 
the solution of systems of linear equations of the form Ri = J2k Pikdk (where 
dk is the layer thickness of the fe-th element phantom in the studied object), 
but with attenuation coefficients pik = p Z^, Pk) that are a priori pre-set 
at several fixed energy values Ei for the basic elements, each with its selected 
atomic number Z^ and density p^. The values of pk = dk/J2k dk depend not 
only on reflexes Ri, but also on the basis chosen. Therefore, such radiogra- 
phy is, in fact, "indirect". Erroneous or inaccurate choice of the basis leads 
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to mistakes (artifacts) in determination of Z e g. For example, a sample with 
low content of a heavy element (low Zh) can be mistakenly identified as a 
sample with a thick layer of a low Zl element. It is clear that, expanding the 
basis and increasing the multiplicity, one can substantially reduce the errors. 
However, even then, errors can be as high as tens of percent. The accuracy 
of traditional methods in 2- and 3-energy radiography allows distinction of 
materials with difference in Z eS by about 50% . So, it is possible to distin- 
guish iron alloys (Z « 26) from wood (Z f» 6) and even from light aluminum 
alloys [Z « 13), but it is not possible to reliably discern the difference be- 
tween water (Z rj 8) and organic materials {Z ~ 6 — 7), iron and many other 
metals, calcium minerals {Z ~ 17) and sand (Z « 12), etc. To distinguish 
among these materials requires accuracy that should be higher by an order 
of magnitude (i.e., not worse than 95% ). 

For a chemical compound with the known formula of the form [X] = 

N 



J' 



.X J 



(here rij is the number of atoms of the i-th kind 

=i 

with atomic number Zj and atomic mass Ai in the complex molecule [X]; N 
is the full number of simple components), Z e g and A e g are determined from 
the expressions 
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where p is an index depending upon the chosen absorption mechanism of 
gamma-quanta. The expression for determination of Z e s is derived from 
the balance of the absorbed energy per each atomic electron of the sub- 
stance. Absorption over all possible physical channels and by all simple 
components of the complex compound is accounted for. Parameter p is 
related to the characteristic dependence of the electron coefficient of lin- 
ear absorption, /i e = fi m /Z = p/pZ (here \i m is the mass attenuation co- 
efficient), on the atomic numbers of simple constituents of the substance 
studied. The choice of p values corresponds to the predominant absorption 
channels at specified radiation energies, i.e. /i e oc Z^ s . In the photo effect 
theory, (/Xe) p hoto °^ Z 3 . In the pair formation theory (/i e )p airs oc Z. For the 
Compton scattering weakly depending upon atomic properties of the mate- 
rial (/i m ) compt oc Z° (/ie) compt Z~ x . In the experiments, such absorption 
character is also well confirmed. Besides this, the parameter p can be consid- 
ered as a fitting parameter giving the best approximation of the absorption 
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cross-section in a specified energy range and for a specified class of sub- 
stances. E.g., in the photo effect energy range intermediate values p = 2 4- 4 
are used. We use for the photo effect p = 3, while for the pair formation 
effect, p — 1, i.e. 



J photo 



N 



N 



riiAiZf / ^ ntAiZi 



.t=i 



8=1 



1/3 



y pairs 



N IN 

^ riiAiZf I n i A i z i 

i=l I i=l 



(2) 

For the Compton scattering, which normally accompanies one of the for- 
mer mechanisms, it should be assumed by definition that p — — 1 and 



N IN 
■^compt = 2^ n iAiZi I fliAi 
i=l I i=l 



(3) 



as the Compton effect cross-section does not depend upon properties of the 
material, but only upon its average electron density. The molar mass Mx = 
J2i n iAi- In practice, relative concentrations q of the simple components of 
the material are often also known. Then, 

n i (where < q < 1), 

instead of (JU)-© we obtain 
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N 



1/p 



A 



eff 



N 
i=l 



(4) 



a=i I i=i 

In Table 1, data are presented on the effective atomic number of selected 
substances that comprise many materials in practical use. The data were cal- 
culated using formulas (JU)-(|I1). It should be noted that the effective atomic 
number is dependent upon the energy range of the ionizing radiation used. 
Its value corresponds to the predominant absorption channel. Using elemen- 
tary inequalities from expressions one can easily obtain the general 
relationship valid for chemical compounds: 



^compt ^ ^pairs — ^photo ) (5) 

which is substantiated by the data of Table 1. For homogeneous mixtures 
(solid, liquid, gaseous), including alloys, these inequalities © can be violated 
at certain concentration ratios of the mixture components. This feature 
(inversion) can be used, for example, in identification of substitution alloys 
or composites. 
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3. Theory and physical model of two-radiography 



Let us consider a simple basic model of two-energy radiography used for 
direct qualitative determination (monitoring) of the effective atomic number 
of a material. A general scheme of such radiography is presented in Fig. 1. 
X-ray and gamma-radiation is attenuated exponentially with linear coeffi- 
cient n = p(E;p; [X]). The latter depends upon the radiation energy E, 
the density of the material p, and its chemical (atomic) composition. Z e g 
is considered as a direct characteristic of the atomic composition. For sim- 
plicity, we assume that: 1) the radiation is monochromatic at the two fixed 
energies E\ and E%\ 2) its spectrum is not changed when passing through 
the object, and 3) scattered radiation can be neglected. By the appropriate 
choice of energy filters, or by using radioactive isotopes as radiation sources 
and subsequent collimation of the radiation in a sufficiently narrow beam, 
rather good approximation of these conditions can be realized in the exper- 
iment. Corrections for these factors can be also accounted for theoretically. 
The non-monochromatic character of the emitted and detected radiation, as 
well as its accumulation due to scattering inside the object (and/or detector), 
generally lower the monitoring efficiency and relative sensitivity. 

Let us write down the (digitized) signal as recorded by detectors in the 
form 



where Vq is the background signal (without an object); dis thickness of local 
cross-section of the object in the direction of ray propagation. One should 
note that the value of Va = Vq (Ei) depends upon full conversion efficiency 
of the system (i.e., the ratio of the energy of the useful electron signal to 
the energy of the initial photon) and upon the radiation source power. Let 
us separate the dependence on the effective atomic number Z eS = Z in the 
attenuation coefficient 



where the functions r (E) , a (E) , x (E) define the energy dependence (as- 
sumed to be universal) of the actual absorption cross-sections for the photo 
effect, Compton effect, and pair formation effect, respectively. It should be 
noted that theoretical determination of Z e g for a complex material using for- 
mulas (HI)-© is based just on the absorption structure as given by (J7J). Let 
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us define the reflex as Ri = R(Ei) = In [Voi/Vj] and go from the system of 
equation JBJ) to the following linear equations: 



R l /{pd)=p mi = [a i Z p + (3 l ]Z . (8) 

Here, we have introduced the mass attenuation coefficient p m = \ij p and 
the monitoring constants a, (3, which depend only upon the radiation energy, 
but not on the properties of the tested material. It is taken into account that 
the normal energy ranges of the photo effect (E < 0.6MeV) and the pair 
formation effect (E > 1.2MeV), rule, sufficiently far from each 

other on the energy scale. Hereafter, we will consider just the two-channel 
absorption mechanism, i.e., either of the photo effect/Compton effect type 
(p = 3) or pair generation/ Compton type (p = 1). Monitoring using only 
the inverse Compton scattering corresponds to p — — 1. In all cases, it is the 
relative reflex of 2-radiography that plays the principal role: 

x _ o m _ \n[V (E 1 )/V(E 1 )] 

X - Rl/R2 = ln[V (E 2 )/V(E 2 )} ■ (9) 

As it follows from equations (jHJ, this value is related to the relative at- 
tenuation coefficient, X = p m (Ei)/p m (E 2 ) = p(Ei)/p(E 2 ). Therefore X 
does not depend upon the geometry (thickness) or the density of the material 
nor does it depend on its other physico-chemical properties, except its effec- 
tive atomic number. In the energy range where total absorption coefficient 
varies monotonically, we have < X < 1, if p{E\) < p(E 2 ), or X > 1, if 
p{E l )>p{E 2 ). 

Solving the system of equations (jHJ) with respect to unknown variables Z 
and 6 = (pd) leads, after reconstruction of the monitoring constants, to the 
expression 



Z p 



Zf + X) n ~ (Z p 2 +1 6 2 ) r 2 ] l[{Z x 9,) ri - (Z 2 9 2 ) r 2 ] 



(10) 



where r\ = C\ 2 R 2 — c 22 R\\r 2 = cuR 2 — c 2 iRi, 9j = pj dj is the surface density 
(its dimensionality [6] = g/cm 2 ) of the j-th reference. Formula (fTTH) involves 
only the radiography data R\ j2 and the calibration data for measurements 
with two reference samples with known (Zj; pj) and fixed thickness dj. The 
calibration data are represented by the matrix Cij = R(Ei\ Zj,8j), where 
i,j = 1,2. In calibration, one should account for the solubility conditions 
of system (jHJ), which are, in fact, conditions on the monitoring feasibility. 
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Hence, it follows that the choice of the radiation range and the references 
should comply with the requirements detc^ = Cnc 2 2 — c 12 c 2 i ^ 0; Z\ ^ Z 2 . 
There are no limitations imposed on 8\ j2 , i.e., there is no limiting relationship 
between density and thickness of the samples. 

Let us present equation (fTUj) in the following more convenient form: 

Z cS = [(k 1 X + k 2 )/(hX + k 4 )] 1/p , (11) 

where the new constants k\, k 2 , k 3 , k 4 are related to the old constants from 
(fTU|) by certain relationships. As the formulations appear to be rather 
bulky, we do not write them down explicitly. The dependence (jllj) is a 
fraction-rational, non-linear, though monotonic, function. Expression (jllj) 
will not be changed under a uniform scaling transformation hi, k 2 , k 3 , k 4 — > 
Mk\, Mk 2 , Mk 3 , Mk 4 (M = const). Therefore, only three of the introduced 
constants are independent. They are determined by calibration using two 
samples of given composition, i.e., by the matrix c^-. However, there could be 
another approach to determination of kx, k 2 , k 3 , k 4 . We use, from the begin- 
ning, the functional relationship (jllj) . The fraction-rational function is unam- 
biguously reconstructed from three reference points. Having carried out the 
measurements for three reference materials with pre-set values of Zi, Z 2 , Z 3 
(note that the density and thickness of the samples are completely arbi- 
trary), we obtain three relationships of the form Zj = fj (X,-; k ly k 2 ,k 3 , fc 4 ). 
The constants are then readily determined from these relationships: 

h = Z\Zl [X x - X 2 ) - Z\Z\ {X x - X 3 ) + ZlZl (X 2 - X 3 ) ; 
k 2 = X\X 2 Z 3 [Z\ — Z 2 ) — XiX 3 Z 2 (Z\ — Z 3 ) + X 2 X 3 Z\ [Z 2 — Z 3 ) ; 

h = x, {zi - zi) - x 2 [z\ - zi) + x 3 (zf - ZD ; 



k 4 = x x x 2 (zf - zd - x,x 3 [z\ - zd + x 2 x 3 (z% - zd , (12) 

where Xj = R\j/R 2 j is the relative reflex for radiography of the j-th 
reference (j = 1,2,3) at fixed energies of 2-radiography. Unlike the case of 
Cjj-calibration for expression (fTUj) . three, and not two, reference samples are to 
be used here. This is related to the fact that density and size (dimensions) of 
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the reference samples are not fixed. Obviously, such a calibration procedure 
is more convenient for experiments and practical applications. 

It is convenient to consider parameter p in Eq. (jllj) as fixed. Earlier, it 
was noted that p = 3 for energies in the photo effect range, and p = 1 for the 
pair formation effect. At the same time, p can be considered as one more un- 
determined constant, especially, if it is not known beforehand which absorp- 
tion mechanisms prevail, or if the energies are used at which all these mech- 
anisms are essential. Than the value of p can be determined numerically by 
one of the approximation methods. It can minimize in the best way the set of 
mean square deviations A n (p) = Z% — [ki (p) X n + k 2 (p)]/[&3 ip) X n + k± (p)] 
for a certain basis of calibration measurements for substances with known 
atomic number values Z n (n = 1, 2, . . .). It should be noted that in the same 
way it is possible to substantially increase the accuracy of effective atomic 
number determination in a class of objects with close Z e s values, for instance, 
to increase the accuracy of distinction between organic substances, etc. 

4. Discussion and experiment data 

The proposed method of "direct" reconstruction of Z e g is free from the 
above-discussed disadvantages related to the monitoring being dependent 
upon the chosen basis. The "direct" method ensures up to 95% accuracy 
for Z e ff determination. In fact, using expressions ©-(^2) to determine the 
relative sensitivity Sz , we obtain an estimate 

S z = (AZ/Z) oc 2 (Ad/d) , (13) 

where AZ is the minimum detectable change in the effective atomic num- 
ber, and Ad are the smallest detectable variations of the object thickness. 
Sensitivity with respect to defect detection Sd = (Ad/d) for multi-energy 
radiography corresponds to the sensitivity of conventional radiography (sep- 
arately for each of the assembly detectors) and is normally of the order of 
several percent. The numerical factor 2 in the expression (|13|) is related to 
the two-energy nature of the monitoring (i.e., sensitivity is assumed to be 
equal for all detectors, Sd,i = Sap)- Consequently, Sz oc 5%, if oc 2.5%. 
This requires spatial resolution of 4 pl/mm . 

The relationship (|13|) determines the accuracy of the proposed method. 
It can be seen that 95% reconstruction of the effective atomic number can 
be achieved for an unknown material. Below, this conclusion is confirmed 
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by experimental data. As noted before, the reason for such dramatic in- 
crease in accuracy is related to the direct reconstruction of from radiography 
measurements without using a fitting procedure by choosing the "basic" ma- 
terials. A formal representation ( "replacement" ) of an arbitrary material by 
superposition of two specified materials with known absorption coefficients 
leads to large errors in determination of Z e g. The existing materials are too 
numerous for such rough approximation for absorption of electromagnetic ra- 
diation. On the contrary, the direct effective atomic number using only the 
relative radiographic reflex © for the inspected material allows us to avoid 
this mistake. It is of an order of (AZ) indirect oc Z/M (M is the order of mul- 
tiplicity), i.e., for indirect reconstruction of Z eS in two-energy radiography 
the errors can reach 50% . Many industrial radiographic installations using 
indirect methods have such accuracy, AZ = ±0.5Z e fj. This is quite sufficient 
to discern organics from inorganics, e.g., to see the difference between heavy 
alloys, Z e s = 26 ± 13, and plastics Z c s = 6 ± 3, but is clearly not enough 
for more accurate measurements in many applications. The direct method 
of Z e ff reconstruction opens here many new possibilities. 

In general, unlike the "synthesis" of Z e g, the direct method is based upon 
"analysis" of the atomic composition. Moreover, in the new approach it is 
possible to limit oneself to the use of just the 2-radiography. This is im- 
portant, because passing over to radiography with higher multiplicity is a 
technically difficult task. To verify the theory, we compared the obtained 
theoretical dependence (jll|) - (|12j) with known experimental data on gamma- 
radiation absorption in a large range of various materials, starting from car- 
bon [Z = 6) and ending with uranium [Z = 92). These results are shown 
in Fig. 2-4. In constructing the theoretical curves, three points were chosen 
as reference ones (defining a fraction-rational function), which corresponded 
to materials with large, small and intermediate Z e ^ values. Three character- 
istic regions of the energy spectrum were considered - i.e., low, middle and 
high energies, corresponding, respectively, to the photo effect, the Compton 
effect and the pair production effect. The data presented are in a very good 
agreement with theory. 

The choice of absorption mechanisms depends upon the radiation energies 
used. The gamma-quanta energy range is mainly determined by the exper- 
iment conditions of character of applications. E.g., in security and customs 
inspection of luggage, medium energies from several tens to several hundreds 
keV are used. This corresponds to the region of combined effects of photo 
absorption and Compton scattering of gamma-quanta in the substance (Figs. 
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2-3). Discarding any of these mechanisms in the medium energy range can 
lead to substantial errors in determination of the effective atomic number. 
Thus, in Fig. 3 a theoretical approximation neglecting the photo effect gives 
much worse results. For inspection of large objects, such as trucks or con- 
tainers, accelerator energies of several MeV are needed. In this case, for 
determination of the effective atomic number one should use a model where 
the pair formation effect and Compton effect are predominant It should be 
stressed that the proposed method can be used with any choice of the ra- 
diation energies. For the data presented in Figs. 2-4, we used only several 
fixed radiation energy values. We have checked that with any other choice 
of these energies the principal law (jllj) remains valid. For the experimental 
data used, deviations from the theoretical dependence did not exceed 5-10% 
, which corresponds to 90-95% accuracy in Z e g determination. 

Detailed analysis of the optimal design for 2-radiography, completion of 
direct experiments, reconstruction of the atomic composition of different ob- 
jects and complex materials, micro- and macro-structures, etc. are the sub- 
jects of our further studies. The approach presented in this paper allows the 
most complete extraction of the information on physical properties of the 
material studied by multi-energy radiography. The number of parameters 
reconstructed corresponds to the radiography multiplicity. (This is a general 
principle of multi-energeticity). This direct analytical approach uses only ab- 
sorption data, showing up the amplitude contrast induced by ionizing radi- 
ation of photons (as monochromatic electromagnetic waves). Together with 
new X-ray analysis methods, including those using synchrotron (coherent) 
radiation exposing the phase contrast, as well as multi-energy holography 
[IB], [H], EH]; micro-focusing and computer tomography [21], [22], etc., the 
proposed approach to monitoring of materials structure widely broadens the 
field of NDT possibilities. This direction is similar to new directions in stud- 
ies of semiconductor and other films. Also, certain prospects are opened for 
application of NDT in development of various microstructures, as well as in 
studies of the distribution profile of implanted nano-clusters of alien atoms 
in the "host" crystal lattice. 

5. Summary 

Thus, a direct approach has been proposed to reconstruction of the atomic 
structure of materials by means of multi-energy radiography. The general ex- 
pressions fllO|) - (fT2"|) obtained for the effective atomic number in 2-monitoring, 
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are closed and convenient for development of computer data processing al- 
gorithms. These formulas can be also used for reconstruction of the struc- 
ture of mixtures, composites, multi-component objects, micro- and macro- 
structures, systems with variable chemical composition, etc. The validity of 
a universal law (fraction-rational dependence with calibration over three ref- 
erence points) in the atomic number radiography is experimentally confirmed 
in different energy ranges of X-ray and gamma-radiation. An essential feature 
of multi-energy radiography is direct extraction of additional information on 
physico-chemical structure of the object under study. This opens new possi- 
bilities in materials science, non-destructive testing and other applied fields 
of science and technology. 
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Table 1. Effective atomic number of various substances with 
respect to the photo effect, zT photo , pair formation effect, zT pairs , and 
Compton scattering, Z compt . 
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Figure captions 

FIG. 1. General scheme of two-energy radiography with reconstruction 
of the effective atomic number of the material. Synthesis consists in mix- 
ing of the fitting basic elements with L - "light", M - "middle" and 
H - "heavy" atomic mass. Analysis is unambiguous reconstruction of 
Z e ff. "Black-and-white" synthesis corresponds to the two-energy radiog- 
raphy, "three-color" scheme (R - "red", G - green, and B - "blue") cor- 
responds to 3-radiography, etc. For the direct method (analysis) proposed in 
this work, 2-radiography is sufficient. 

FIG. 2. Dependence of the effective atomic number upon 2-radiography 
reflexes in the regions of photo effect and Compton scattering. Theoreti- 
cal dependence (solid line) and experimental points for materials of known 
composition are indicated. Geometry (dimensions) of samples is arbitrary. 

FIG. 3. Dependence of the effective atomic number upon 2-radiography 
reflexes in the region of intermediate ionizing radiation energies. Theoretical 
curves are presented modeling the predominant role of one of two possible 
absorption channels. The best agreement with experimental data is obtained 
when a mixed absorption mechanism (photo effect /Compton effect) is chosen. 

FIG. 4. Dependence of the effective atomic number upon 2-radiography 
reflexes in the region of Compton scattering and the pair formation effect. 
When the detected energy ranges are moved apart, sensitivity of the method 
is not worsened. 
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Figure 3: 
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Figure 4: 
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